Correlated ground states with (spontaneously) broken time-reversal symmetry 
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We propose a self-consistent scheme for the determination of the ground-state (GS) properties of 
interacting electrons in a magnetic field, and of systems whose GS's time-reversal-symmetry (TRS) 
is spontaneously broken. It is based on a newly-developed many-body perturbation theory that is 
valid, irrespective of the strength of correlation, provided the GS number densities »i|(r), nj^(r), 
and the totaJ paramagnetic particle flux density are pure-state non-interacting u-representable. Our 
approach can in particular be applied to (modulated) two-dimensional electron systems in the frac- 
tional quantum-Hall regime. 
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Consider the following Hamiltonian that in the non- 
relativistic limit governs the behaviour of a general sys- 
tem of interacting electrons (we employ the SI units): 

H^Y. J d'r (r) 1^ [-zW + eA{r)f 



+ v{v) + -gtiBB{v)a,{a)'^r^A^) + V, (1) 

with V={\I2)Y.,^^, j d^rd^r' i^l{v)4,l{v') v,{r ~ r') 

V'(j'(r')'0cr(r) the interaction part (the quantities in these 
expressions are defined in Ref. [|l| ) . Here for simplicity we 
have assumed that the magnetic field B{r) = V A A{r), 
with A the external vector potential, points in a definite 
direction, i.e. B{r) = B{r)ez- 

Efforts for determining the electron-electron interac- 
tion effects on the physical properties of systems de- 
scribed by H include, e.g., use of the many-body 
perturbation theory at some level of approximation, 
various quantum Monte-Carlo techniques and exact- 
diagonalisation methods. Application of the (ordinary) 
perturbation theory is considered to be limited to so- 
called 'weakly-correlated' systems, and although suitable 
for dealing with 'strongly-correlated' systems, the lat- 
ter two categories of mentioned methods are applica- 
ble only to systems consisting of relatively small num- 
ber of particles or those whose Hilbert space is small. 
In this Communication we introduce a self-consistent 
many-body perturbation theory that, in principle, can 
be applied to weakly as well as strongly correlated sys- 
tems. Here our attention will be mainly focussed towards 
some GS properties. In order to be able to describe GS's 
with spontaneously broken TRS, we retain A in Eq. (|l|) 
and identify the case of zero external magnetic field with 
B{r) 0. 

In a previous work we have discussed a fundamental 
problem from which any perturbation theory can suffer: 
that despite the possible convergence of a perturbation 
series, the ultimate results may not even approximately 
be related to the quantities of interest. The reason for 



this type of breakdown of the perturbation theory lies in 
that the GS of the non-interaction Hamiltonian Hq=H — 
V (possibly modified by some effective one-body term), 
may not be adiabatically connected with that of the fully 
interacting system Q . 

Our analysis in Ref. shows that a perturbation the- 
ory based on a non-interacting Hamiltonian Hq whose 
ground state \%) satisfies {%\di\%) = (^"010,1*0) 
=0^, z 1, 2, . . . , A/", with l^-o) the GS of the fully inter- 
acting system and {Oi\i = 1, 2, . . . ,N'}=Sj\/ a specified 
set of quantities that uniquely determine the ground state 
of the system (see following paragraph), is an uncondi- 
tionally valid perturbation theory [Q. We observe that 
for Oi — (^ojOil^Q) to be satisfied it is necessary that Oi 
be pure-state non-interacting w-representable (for defini- 
tion see Ref. 

For the work that we present in this Communication 
we rely on a theorem due to Vignale and Rasolt (VR) 
Q according to which the GS of the Hamiltonian given 
in Eq. (^ is a unique functional of the number densities 
ncr(r), a S {Tji}; and the total paramagnetic particle 
fluxdensity jp(r)=_7p.|(r)-t-jp.|(r). If the external vector 
potential A were spin dependent, denoted by A^, then 
the GS would become a unique functional of (r) and 
j„{r) for a i- To keep our approach general, so that 
both and Jp.|, pertaining to the interacting system, 
can be calculated ^ , in what follows we formally assume 
to have a spin-dependent external vector potential; for 
the actual calculations we take A-^ = A| = A. Thus for 
the present case we have = 84 = {no-, jp.a-\ c —1, 1} 
in which n^(r) =iijl{r)%jj^{r), and jp.„{r) ^{-in/[2me]) 
{i/^t(r) [VV^^(r)] -[VV^J(r)] ^^<x(r)}. The Hamihonian 
Hq for the case at hand coincides with the Kohn-Sham 
(KS) Hamiltonian as introduced by VR ||]. We have 

HkS = EaJd^'r i>l(v)HKS;a{v)^'r{v), in which § 
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^ (-zW + e[A(r) + A,e;a(r)])' 
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+v{r) + -gfiBB{r)<72{a) + ^^(r) + v^c-.^ir). (2) 

Let now Ga{j^t,r't') denote the single-particle Green 
function corresponding to H, and Gks-.u that corre- 
sponding to Hks- Making use of the results 



njr) - 



n^(r) = -iG^{Yt,Tt + Q+), 

jp-A^) = li: l™r'^r (V - V) G,(rt, v't + 0+), 



(3) 



which, provided n„ and j^.^ are pure-state non- 
interacting ti-representable, are by construction also valid 
for Go- Gks:<7, and a perturbation expansion for G^ in 
terms of Gks-.o-, we obtain for a c?-dimensional system a 
set of 2(d+l) coupled non-linear equations which we have 
diagrammatically represented in Fig. 1. From these equa- 
tions the two {d+ l)-vectors {v^c-a, A^c-a)^ <^ ^ {T,i}, 
can be determined and consequently the Hks on the ba- 
sis of which an unconditionally valid perturbation theory 
can be set up The unfamiliar Feynman diagrams in 
Fig. 1 have their origin in the specific choice for the "un- 
perturbed" Hamiltonian. For Hi = H — Hks we have 

-Y, I d\ (r) {^;^(r) -f v,,,„{v)} V^<,(r) + V. (4) 

For casting Hi into a form that makes application of the 
standard procedures of the many-body perturbation the- 
ory [|[ possible, we define the following operator-valued 
non-local potential: 

1 pn 

u^^.^r, r')=— <5(r - r')A,,;,(r) • (V - V) . (5) 

Through this the first term on the right-hand side of 
Eq. (U) can be written as [— X^o- / <Pt(Pt' Ua:c;CT(r, r') 
V'J(r)i/)cr(r')]- In Fig- 1 the directed wiggly hues, pointing 
from r'fj to rcr, stand for {—i/Ti)[—Uxc;<j{'^i'^')]- For the 
rules concerning evaluation of contribution of diagrams 
see Ref. [|. 

Earlier Sham [||, by imposing G — Gks at r' — r, 
t' = t + 0^ for systems of spin-less electrons (with no 
broken TRS), obtained an implicit expression for Vxc- 
Sham, however, did not address the problem concerning 
the validity of the many-body perturbation theory which 
is central to our present work. 

Let the contributions of diagrams (ai), {hi) and (ci) 
in Fig. 1 be denoted by n^a\Y) and j^^\{v), with x = 
oi , 6i , ci . The non-linear equations associated with these 
first-order diagrams read (with cr G {t, i}) 
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FIG. 1. Diagrammatic representation of Eq. (p|) with Go- 
(thick line with single arrow) and GkS;ct (double-arrowed thin 
line) on the right-hand sides (upper parts of (A) and (B)). 
Contributions of diagrams (ai), (bi), (ci), ... add to zero 
for the correct v^c-.a (wiggly line with one loose end; the as- 
sociated expression equals (— j/?i)[— fxcio-]) and Axc;cr which 
defines Uxc;tT (directed wiggly line) according to Eq. (|^). The 
pulse-train-like line in (B) implies the non-local operation 
{-h/[2m^]) hmr'_r (V - V) in Eq. (||). Lower parts of (A) 
and (B) correspond, respectively, to the first and the second 
equation in Eq. Note that in first order only one diagram 
(i.e. (ci)) explicitly depends on the Coulomb interaction, rep- 
resented by the broken line (corresponding to {—i/h)vc); this 
line can be taken to represent also the dynamically screened 
interaction function. 

ni"^)(r)+n^)(r)+n^)(r) = 0, 
4':^Hr)+4^i)(r)+jti)(r)=0, 

in which 

n^f\r) ^ -2Rc Ef> uj^^j^, g.,, ,(r)} , 
j(^(r) ^ -2Re C^S.^'^'^^^'.^ ^ 

Ui.:l=ld'r' f.c;a(r')e.;.,.'(r')/(e.;.' 

Ui.ll=ejd^r' A,e;,(r') • J,;,,^ (r')/(ea;s' - e.;.), 

X [Jd^r" v,{v' ~ v")Q,,s".,Ar")\ /(e,;,, - e,;,), 

6'<T;s,s'(r)=C;s(l")V'ff;s'(r); J'^Tis^s' (r) = (-^V [2"^e] ) 

x{V':;,(r)[VV'.;.'(r)] - [VC;s(r)]^a;.'(r)}. 

(8) 

We have n„(r) = YA Q'y.sA'^)^ JpA'^) = J^-sA^)- 
For A = Aj:c = 0, the first equation in Eq. (|^) reduces 
to one embodying the 'optimised-potential method' [Eo|. 



(7) 
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Above (E. ) stands for ErJi iET=N.+i) and 
4>a;s denotes an eigenfunction of HkS;<j in Eq. (g), 
with ea;s the corresponding eigenvalue {ea;s < ecr;s+i)- 
Thus, Eq. not only explicitly depends on Vxc;a and 
Axc-a (via Eqs. (^), but also impHcitly, through 
{^/'cr;s,£(j;s}. Since in solving Eq. (^) we exactly di- 
agonalise the KS Hamiltonian corresponding to any 
{vxc;cr,-A.xc;a), the self-consistcnt {vxc;a,Axc;a) takes ac- 
count of the electron-electron interaction effects, in so 
far as present in the employed diagrammatic expansion 
for Go-, to infinite order. Moreover, it can be shown 
(cf. Ref. [^) that a solution {vxc-a, -^xc-a) of Eq. (||) 
based on a finite-order perturbation expansion for 
annihilates the combined contributions to Un and j^.^ of 
all higher-order diagrams that have a common part in 
addition to one of the lower-order number-density dia- 
grams (Fig. 1 (A)) that have been taken into account. 
Note that in solving Eq. (||), and N^^ satisfying 
N-\ -\- Ni = Ne, must be calculated self-consistently: for 
a given N^., 7V|, or N^, is determined by the require- 
ment that £|;i, . . . ,ef;Ar^, ei-i, . . . ,si-Ni are the lowest 
Ne eigenvalues of {HKS;a}- 

The non-linearity of Eq. (^ implies that there is a mul- 
tiplicity of solutions. The uniqueness theorem of VR ||] 
implies, however, that only one of the solutions corre- 
sponds to the GS (solutions a + Vxc;'\, ct + Vxc;i, with a a 
constant, are identical). It is important that all solutions 
of Eq. (^) correspond to some ||] eigenstate of H. To sin- 
gle out the GS solution {vxc;a , Axc-a) , let the GS of the 
system corresponding to the external potential v — vq 
and the external magnetic field strength B = Bq (and 
therefore Aq, in some gauge) be known. For the actual 
potential v and magnetic field strength B we define: v^'^'^ 
=vq +X{v - Wo) and A^^'^ =Ao +X'{A - Aq). For de- 
termining the GS solution of the non-linear equations we 
choose some trajectory on the (A, A')-plane connecting 
(0,0) with (1, 1). By starting from (0,0), at small steps 
along the trajectory we solve Eq. (^) to self-consistency. 
In doing so, at each step we take the SC solution of 
the immediately earlier step as the trial solution. In the 
event of encountering degeneracy (level crossing at, say, 
(Ao, A'o)), we can easily select out the GS solution by re- 
alizing that at the point of degeneracy the energies of the 
degenerate states have discontinuous derivatives and the 
jump in the derivative of the GS total-energy curve is al- 
ways negative and has the largest magnitude amongst the 
jump values corresponding to other energy curves. Now 
since the solutions found pertain to eigenstates of the 
Hamiltonian corresponding to v^^^^°^ and A*-^ °\ we 
can apply the Hellmann-Feynman theorem 1^, |^ which 
provides us with the derivatives of the eigenenergies with 
respect to A and A'. A simple calculation shows that 
ni^\ n\^^ and jl^J + j'i^i'' are sufficient for determining 
the mentioned derivatives Note that through inte- 
gration of the available total-energy derivatives along the 



chosen path, the GS total energy is readily obtained. 




FIG. 2. Particle number and flux densities corresponding 
to a two-electron quantum dot in the (M^, Sz) state with 
and Sz, respectively, the orbital- and spin-angular momentum 
quantum numbers along the quantisation axis. The param- 
eters chosen are those employed in Refs. |Q, B = 1.0 
T, v{r) = with m* = 0.067me and = 3.37 meV; 
e,- — 12.4, g — —0.44; see Ref. jl|; the effective Bohr radius ao 
equals 9.794 x 10~^ m. The current density —ejp;a{r), with 
a =1, is in units of 2.95 x 10^ A/m (note the — e); jp;a, the 
azimuthal component of jp-.a (the radial component is vanish- 
ing), corresponds to the symmetric gauge, A = B A r/2. 

Our proposed framework yields and jp.^ that ex- 
actly satisfy the static equation of continuity. This is 
because both of these are derived from the same single- 
particle Green function, Gks;^- Further, provided the 
'associated' diagrams (see further on) corresponding to a 
particular order of the perturbation expansion for Ga are 
taken into account, our framework is also gauge invariant. 
To see this clearly, let the 4-potential {vxc;a, Axc-a) cor- 
respond to (w, A^). For definiteness suppose we have ob- 
tained the former by taking into account diagrams (ai), 
(6i) and (ci) in Fig. 1, or by solving Eq. (||). Since Aa- 
A „ = A„ -t- VAo- in a gauge-invariant theory must lead 
to n„ and j^.^ j^.^ - {e / me)naV ■, it follows 

that for v'xc:a and A'^^.^, corresponding to A' cr, hold [||: 
'v'xc-a = '"xc-a +{e^lm^) Axc-a 'V A^- and A^^.^ = Axc-^. 
Our framework would be gauge non-invariant if Eq. (|g) 
would not be satisfied by {v'xc-cy^ -^xc-.a)- We now show 
that this is not the case. First, eo-;s does not depend on 
the choice of gauge. Further, from the explicit expres- 
sions in Eq. (|^) it can be shown that Qa-^a,s'^ similar to 
no-, is gauge invariant, and that for A„ — > A o-, J a-.s.s' 
-> Jcr;s,s' -{e/me)Qa-s,s' VA^-. It foUows from Eq. (|) 
that u'^l^, +ul''.l\, and ul'^l^^, are gauge invariant (dia- 
grams (fli) and (6i) in Fig. 1 are 'associated'). Thus n'^^^ 
+71^^'' and ni'^^' are gauge invariant and j^^^^ +jp^ ~^ 
+3f-^) -(e/me) [n^^^^ W^'^) V A,, ^'r^) j"^ 
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— (e/TOe) n-a V Ao-. Therefore satisfaction of Eq. (|^) 
does not depend on the choice of gauge. 

We have apphed our formahsm to a cyhndrically sym- 
metric quantum dot with a paraboUc confining potential, 
taking into account only the first-order diagrams that are 
explicitly shown in Fig. 1. In Fig. 2 we present the calcu- 
lated electronic number and flux densities in the GS's of 
definite symmetries (see caption) and compare the for- 
mer with some available results [|l^, — note that 
= in the (M^ = 0,S^ = 0) GS. It turns out that in 
the (1, 1) state, our calculated n [= +ni; for the (1, 1) 
state, n = n|] is, within the numerical accuracy of the 
calculations, identical with the exact n. On the other 
hand, for small values of r our (0, 0) GS n(r) overesti- 
mates the exact n{r) by several percents. Both of these 
density profiles are almost identical with the Hartree- 
Fock (HF) results ||l2|. There is however a fundamental 
difference between the HF results and those according to 
the present scheme. According to the HF approach, for 
all B ^ 0, the GS of the system under consideration is 
a (1, 1) state, in obvious contradiction with the exact re- 
sults For instance, for B = IT, the exact GS is a 
(0,0) state, and the (first) transition to the (1,1) state 
takes place at B = 2.0T. In agreement with this, we find 
for _B = IT the lowest-lying state to be a (0, 0) state. It 
should be mentioned that the results in Fig. 2 labelled 
by (1, 1) correspond to the lowest-energy (1, 1) state; for 
B = IT, as indicated, the state corresponding to the ab- 
solute minimum of energy is a (0, 0) state. It is therefore 
important to emphasise that the VR theorem |^ is also 
valid for excited states which are minimum-energy states 
corresponding to specified symmetries. This follows from 
the fact that the variational principle, which underlies the 
VR theorem, can also be applied in symmetry-restricted 
Hilbert spaces. For comparison, in Fig. 2 we present the 
results obtained within the local-density-approximation 
(LDA) scheme |l|l. 

Concerning the overestimation in the vicinity of r = 
of our calculated n{r) corresponding to the (0,0) state, 
our preliminary calculations indicate that this is substan- 
tially suppressed through replacing Vc in diagrams (cl) of 
Fig. 1 by the dynamically-screened interaction function 
within the random-phase approximation. 

We have analysed the asymptotic behaviour of the 
fimctions Vxc:a{i^), ^xc;CT(r), for |r| — > 0, oo, correspond- 
ing to the system under consideration, both within the 
framework of the present formalism and that of the LDA. 
The results will be reported elsewhere jl^]. We only 
mention that unless appropriate measures are taken, the 
current-carrying GS's of the LDA |^ are unstable. 

In conclusion, we have introduced a self-consistent per- 
turbation theory for interacting electrons in presence or 
absence of an external magnetic field. In the latter case 
the system can possibly have a spontaneously broken 
TRS, such as is the case in open-shell atoms. Already 



the first stage in the application of this theory provides 
us with the scalar and vector exchange-correlation poten- 
tials that determine the Kohn-Sham Hamiltonian within 
the framework of the current-and-spin density-functional 
theory. This Hamiltonian forms the basis for construc- 
tion of reliable perturbation expansions for various quan- 
tities, including those corresponding to the excited states 
of the interacting system, such as energies of the ele- 
mentary excitations. We propose use of our method for 
determining properties of (modulated) two-dimensional 
electron systems in the fractional quantum-Hall regime 
where electrons are strongly correlated. Work in this di- 
rection is in progress. 
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